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We present a correlative study of structural and optical properties of natural defects in planar 
semiconductor microcavities grown by molecular beam epitaxy, which are showing a localized polari- 
ton spectrum as reported in Zajac et al, Phys. Rev. B 85, 165309 (2012). The three-dimensional 
spatial structure of the defects was studied using combined focussed ion beam (FIB) and scanning 
electron microscopy (SEM). We find that the defects originate from a local increase of a GaAs layer 
thickness. Modulation heights of up to 140 nm for oval defects and 90 nm for round defects are 
found, while the lateral extension is about 2/xm for oval and 4/im for round defects. The GaAs 
thickness increase is attributed to Ga droplets deposited during growth due to Ga cell spitting. Fol- 
lowing the droplet deposition, the thickness modulation expands laterally while reducing its height, 
yielding oval to round mounds of the interfaces and the surface. With increasing growth tempera- 
ture, the ellipticity of the mounds is decreasing and their size is increasing. This suggests that the 
expansion is related to the surface mobility of Ga, which with increasing temperature is increasing 
and reducing its anisotropy between the [110] and [110] crystallographic directions. Comprehensive 
data consisting of surface profiles of defects measured using differential interference contrast (DIC) 
microscopy, volume information obtained using FIB/SEM, and characterization of the resulting 
confined polariton spectrum are presented. 



I. INTRODUCTION 

Fundamental physics was demonstrated in planar semi- 
conductor microcavities over last 10 years including 
Bose-Einstein condensation of exciton-polaritons 1 , for- 
mation of vortices 2 and superfluidity^. To understand 
these two-dimensional inhomogeneous non-equilibrium 
systems, it is important to understand and harness spa- 
tial disorder in these structures. A significant con- 
tribution to polariton disorder in molecular beam epi- 
taxy (MBE) grown microcavities is photonic disorder 
including a cross-hatched dislocation pattern^ 5 and 
point-like-defects (PDs) 6 . Such disorder creates a po- 
tential landscape for the in-plane motion of polaritons, 
which results in inhomogeneous broadening, enhanced 
backscattering 7 8 , and creation of localized polariton con- 
densates and polariton vortices^. 

In this work, we present a correlative study of the 
structural and optical properties of natural defects in pla- 
nar semiconductor GaAs/AlAs microcavities grown by 
MBE. The paper is organized as follows: in Sec.[ll]we re- 
view the literature on defects in MBE grown GaAs struc- 



tures, in Sec. Ill we discuss the samples and experimental 
methods used, in Sec. |IV| we present the experimental re- 
sults, followed by a discussion of the defect formation in 
Sec. M and conclusions in Sec.lVll 



II. 



POINT-LIKE DEFECTS IN GAAS 
HETEROSTRUCTURES 



Point-like-defects observed in MBE-grown GaAs based 
heterostructures are classified into several types, with 
the most common being oval defects. It was reported 
that oval defects originate from an ex cess of Ga, Ga 
droplets, or surfa ce contamination^^, or low growth 
temperature d 12 * 13 l Oval defects were extensively stud- 



ied in the 1990s as they were responsible for the failure 
of electronic devices such as field-effect transistors.^ The 
defects have typical sizes in the order of several /mi and 
a roughly 3:1 aspect ratio along the [ll0]:[110] crystal 
directions.^ Their height on the surface were found to 
be several tens of nanometers. Another type of defect 
observed in this work were round defects having simi- 
lar diameters and heights as oval defects. They were 
attributed 14 to Ga oxide or effusion cell spitting. We 
found that the defects investigated in our work originate 
from a GaAs thickness modulation, which we attribute to 
Ga droplets with sizes in the order of 100 nm emitted by 
the Ga source during growth. The droplet formation was 
previously ascribed 10 to an inhomogeneous temperature 
distribution in the Ga crucible of the effusion cell. Specif- 
ically, Ga cools near the orifice of the crucible and, since it 
does not wet the pyrolytic Boron Nitride (PBN) crucible 
surface, forms droplets which can fall back into the liq- 
uid Ga, causing a spatter of smaller Ga droplets. Recom- 
mended methods to reduce this mechanism include 1) use 
solid instead of liquid Ga; 2) modification of the orifice 
geometry of the Ga cell to inhibit condensed Ga droplets 
entering into the Ga source; 3) creating a positive axial 
temperature gradient toward the orifice to prevent con- 
densation of Ga; 4) treating the crucible with Al, which 
forms a A1N layer which Ga is wetting, suppressing the 
formation of droplets. Another mechanism for formation 
of oval defects is suggested in Ref.[T5l and referred to as 
Ga source spitting. During heating of the Ga source up to 
1200°C, within the range of Ga evaporation^, explosions 
in the Ga liquid were observed, which resulted in Gal- 
lium droplets deposition on the MBE chamber walls. It 
was speculated that these explosions were due to Ga203 
shells encapsulating Ga and creating an effusion barrier. 
Another possible mechanism 1 ^ is that particulates re- 
leased from the walls of the MBE chamber are entering 
the molten Ga in the crucible causing a turbulent re- 
action. Summarizing, a number of mechanisms for the 
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Sample 
cavity length 

DBR periods top(bottom) 



growth 

temperature 

(°C) 



DBR AlAs 
DBR GaAs 
cavity GaAs 



MCI 
lAc 
24(27) 



MC2 
2A C 
23(26) 



715 
660 
630 



590 
590 
590 



TABLE I: Parameters of samples MCI and MC2. 




FIG. 1: Sketch of the optical imaging spectroscopy setup 
used to measure the localized polariton states. Ml: Gim- 
bal mounted mirror, L1-L5: Lenses, MC: Microcavity sample, 
LS1,LS2: movable lenses for imaging, dashed lines: removable 
mirrors, BS: Beam-splitter. 



Ga nano-droplet formation have been suggested, and the 
mechanism dominant in a given growth is not obvious. 



III. SAMPLES AND EXPERIMENTAL 
METHODS 

In this work we investigated two microcavity samples, 
MCI and MC2, grown in a VG Semicon V90 MBE ma- 
chine with a hot-lip Veeco 'SUMO' cell as Ga source, with 
structures given in Table [T[ Sample MCI was studied in 
Ref.lHl During the growth of MCI, the wafer temperature 
was ramped up to 715 °C for the AlAs Bragg layers and 
down to 660° C for GaAs Bragg layers, while the cavity 
layer was grown at 630 °C. During the growth of MC2 in- 
stead, the growth temperature was 590 °C for all layers. 
The two samples show a significantly different aspect ra- 
tio of defects on their surface. In MCI they are essentially 
round (see Fig.[4|, while in MC2 they have an aspect ra- 
tio between 3:1 and 2:1 along the [110] : [110] direction as 
shown in Fig.[8j At a temperature of T=80 K, the cavity 
mode energy in the center of the wafer of MCI (MC2) 
is at huj c = 1.480(1.431) eV, respectively, while the bulk 
GaAs exciton resonance of the cavity layer is at 1.508 eV. 

The low temperature optical measurements were per- 
formed using the optical setup sketched in Fig. [I] The 
samples were mounted strain-free on a mechanical trans- 
lation stage moving along the sample surface (x, y) in 
a bath cryostat at T = 80 K in nitrogen gas at 100- 
300mbar. Two aspheric lenses of 8 mm focal length and 
0.5 numerical aperture (NA) were mounted at the oppos- 
ing faces of the sample inside the cryostat to focus the ex- 
citation and collimate the emission, respectively, provid- 
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FIG. 2: SEM images of oval defect PD3. a) prior to milling, 
green lines indicate defect extension on the surface, b) af- 
ter milling the well exposing the epilayer cross-section at its 
side wall. On the lower part of the images, the edge of the 
alignment photomask is visible. 



ing a diffraction limited resolution of 1 fim. The axial po- 
sitions of both lenses were adjustable at low temperatures 
to control the focus of excitation and detection. The ex- 
citation was provided by a mode- locked Ti: Sapphire laser 
(Coherent Mira) delivering 100 fs pulses at 76 MHz repe- 
tition rate and a spectral width of approximately 20 meV. 
The transmission of the samples excited from the sub- 
strate side and detected from the epi-side was imaged 
onto the input slit of an imaging spectrometer with a 
spectral resolution of 15 fieV. Scanning of the sample im- 
age across the spectrometer input slit while keeping the 
directional image on the spectrometer grating fixed for 
two-dimensional hyperspectral imaging was achieved by 
moving the lenses LSI and LS2 appropriately^. 

Spatial height profiles of defects on the sample sur- 
face were measured with differential interference contrast 
(DIC) microscopy using an Olympus BX-50 microscope 
with a 20x 0.5NA objective. DIC images in green light 
(wavelength range 525-565 nm) were taken by a Canon 
EOS 500D camera with an array of 4752 x 3168 pixels 
of 4.8 /im size in the intermediate image plane. The re- 
sulting image resolution was 650 nm in the plane of the 
sample, and about 2nm in the plane perpendicular to 
the sample surface using quantitative DIC. Details on 
the procedure used to extract height profiles of defects 
using DIC microscopy are given in the Appendix. 

To investigate the sample structure below the sur- 
face, we used a Carl Zeiss XB1540 Cross-Beam focussed- 
ion-beam (FIB) microscope which combines an ion- 
beam milling/imaging column with field-emission scan- 
ning electron microscope (FESEM) 19 , providing an imag- 
ing resolution of about 50 nm for the samples studied. 
The internal epi-layer structure was exposed by FIB 
milling. Smooth cross-sections were obtained using a 
two stage milling procedure. In the first step, a high 
beam current of 2 nA, was used, resulting in fast milling 
but leaving a rough and inhomogeneous interface due to 
sputtering and redeposition of material. In the second 
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FIG. 3: Hyperspectral imaging of polariton states bound to PD1 in MCI, measuring the spatially and spectrally resolved 
transmission intensity I(x,y,uj). First and fourth column: l(0,y,u). The energy is shown relative to the polariton band edge 
at hoj c = 1.4789 eV. Intensity on a logarithmic color scale as indicated. Second and fifth column: Real space intensity maps 
of individual states I(x,y,uJi). The state number / and the relative energy h(uoi — uj c ) are given, and the orange lines indicate 
the related peaks in 7(0, y, cj). Third and sixth column: Real space distributions |^ n ,m,g| 2 of corresponding Mathieu functions, 
labeled by their parity ('e': even, 'o': odd), radial (n) and angular (m) order, and parameter q (see Eq.|l]) and Eq.Q). Color 
scale as for measured data indicated in the first column. 



step, the surface was polished with a lower beam cur- 
rent of 200 pA removing a layer of about 500 nm per cut 
resulting in a negligible surface roughness. In the next 
steps, layers with a thickness of 1 jam or less were re- 
moved using the low beam current. Different stages of 
this milling procedure are shown in Fig. [2] Before milling, 
the oval defect (PD3) is seen in SEM (Fig.[2}i), with the 
vertical image scale corrected for the viewing angle of 
36° to the sample surface. A rectangular well of about 
14 jam width and 7 jam depth is milled into the surface to 



one side of the defect (Fig.^), exposing a cross-section 
through the epi-layers to be measured at its side walls. 
After imaging the exposed cross-section with SEM, the 
subsequent cross-section at a controlled distance further 
into the structure is milled. Steps sizes of about 1 jam 
were used at the outskirts of the defect, reducing down 
to 100 nm at its center to resolve the defect source. The 
cross-sections were 6 — 7 jam deep to expose the complete 
epitaxial structure, and their width was adjusted to the 
lateral extension of the defect observed at the surface. In 



order to mark defects on the sample surface for the cor- 
relative studies, a gold alignment mask was fabricated on 
the surface by photolithography. The mask consisted of 
a grid of 400 /im x 400 /im squares with column and row 
indexing. Considering the defect density in the order 
of 10 3 /cm 2 , this mask allows to trace individual defects 
through the different measurement techniques used in the 
present investigation. 



IV. RESULTS 
A. Sample MCI 

The localized polariton states in the round defects of 
sample MCI were examined previously in Ref. 6, Here 
we report on the correlation between the states and the 
three-dimensional structure of these defects using two 
defects referred to as PD1 and PD2 as examples. The 
localized polariton states of PD1 are observed in the hy- 
perspectral transmission images shown in Fig. [3] The 
defect shows 16 localized modes down to — 9meV below 
the polariton band edge. It can be noted that the lowest 
observed stated is of p^-type symmetry. We expect to 
have 2 states with lower energies, another p^-state and 
an s-state. These were not recorded. 

To qualitatively understand the states bound to 
this defect, we compare them with solutions of 
the two-dimensional time-dependent wave equation 
(pi + dy — d 2 /v 2 ) ^ = with the velocity v for elliptical 
boundary conditions. Using elliptical coordinates x = 
ccosh (£) cos (77), y = csinh (£) sin (77) with the focus dis- 
tance c, and the ansatz ^(x,y,t) = R(£)<&(rj)e lU}t results 
in the ordinary and modified Mathieu equations for the 
angular part $(£) and radial part respectively!^ 



+ A- 2g cos (277))$ = 
(<9f + A-2<2cosh(277)) R = 



(1) 
(2) 



where q = (ouc/2v) 2 is the square of the normalized fre- 
quency. The solutions of $ and R are angular and radial 
Mathieu functions. For a given q, the periodic bound- 
ary condition of the angular part $ in Eq.Q deter- 
mines a series of A m (g) with ascending number of nodes 
m = 0, 1, 2, each of which except m = is a doublet, 
having an odd (o) or even (e) symmetry for inversion of 
y. The elliptical boundary is given by a unique c and £0? 
at which the boundary condition, for example i?(£o) = 0, 
holds. This condition and Eq.Q using A = A m (g) deter- 
mines the mode frequencies q njTn corresponding to modes 
with ascending number of nodes n = 0, 1, 2, .. in radial di- 
rection. Analytic expressions for the solutions are given 
in Ref.EH 

Since the polaritons show a quadratic in-plane disper- 
sion for small wavevectors, their in-plane motion is well 
described by the Schrodinger equation.^ By modifying 
the definition of the above solutions for the Helmholtz 
equation are also solving the Schrodinger equation for el- 
liptical boundary conditions, such as the elliptical quan- 
tum well with infinite barriers. 23 

Another family of analytic solutions of the Schrodinger 
equation with elliptical symmetry are given by 
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FIG. 4: Structural characterization of PD1. a) DIC image 
of the sample surface (linear grey scale) and resulting height 
profile across the defect center (green line), b) SEM images 
of cross-sections through the epitaxial structure, taken along 
the black lines indicated in a). The relative distances from 
the edge of the defect are indicated. 



the Her mite- Gaussian modes of an anisotropic two- 
dimensional harmonic oscillator. However, we found 
that they do not describe the observed distributions well, 
since, as we will see later, the confining potential of PD1 
is not similar to a parabolic potential but rather to an 
elliptical well. 

The solutions 2 -^ \^ n ,m,q\ 2 with the mode orders n, m 
assigned to the measured states of intensities \ipi\ 2 are 
shown in Fig.[3J yielding a qualitative agreement. For 
each energy state, q was adjusted in order to reproduce 
the experimental patterns. 

The structural characterization of PD1 by DIC and 
FIB/SEM is given in Fig.[4] The DIC data shows a di- 
ameter of the defect on the surface of about 6 /im, simi- 
lar to the extension of the spatially resolved transmission 
from this defect, and a height up to about 15 nm. The 
FIB/SEM cross-sections show the origin of the defect - 
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FIG. 5: Hyperspectral imaging of polariton states bound to 
PD2 in MCI with huj c = 1.4787 eV, detailed description as for 
Fig.0 



a thickened GaAs layer with a center depression in the 
third DBR period above the cavity, extending over 4 /im 
in and adding about 60 nm in thickness in the center, 
visible in the y = 4.5 /im cross-section. Further discus- 
sion of the defect growth dynamics will be given later. 

We now move to the second defect PD2 on MCI, for 
which the hyperspectral transmission images are shown 
in Fig. [5] The states are arranged along a ring of about 
6 fim diameter, with the lowest state localized at small y, 
and higher states gradually extending along the ring, as 
in a one-dimensional harmonic oscillator, until the whole 
ring is filled. The underlying potential for the polaritons 
appears to be a ring-shaped well, with a depth decreasing 
with increasing y. The analysis of the potential from the 
states shown in Sec. IV C] confirms this interpretation. 

The structural characterization of PD2 by DIC and 
FIB/SEM is given in Fig. [6] The surface profile is similar 
to PD1, with a size of about 7 jam diameter. The defect 
source is much deeper in the structure than for PD1, in 
the 23rd period of the DBR below the cavity, and has 
a larger height of about 90 nm. This height is about 
twice the nominal height of the GaAs A c /4 layer, and 
leads to discontinuities of the DBRs for about 4 layers 
above the defect. The deep center depression of the defect 
is consistent with the ring-shaped polariton confinement 
potential shown in Fig. ITT] 
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FIG. 6: Structural characterization of PD 2, detailed descrip- 
tion as for Fig.[4] 



B. Sample MC2 

This sample was grown at lower temperature than 
MCI (see Table [l|, and shows oval-shaped defects. Two 
examples of defects for this samples are given here, re- 
ferred to as PD3 and PD4. The polariton states of PD3 
are visible in the hyperspectral transmission images in 
Fig. [7] The states come in nearly degenerate pairs with 
point reflected wavefunctions, extended along the x-axis 
([110]), for example the state pairs (1,2) and (3,4), (5,6). 
The states also have an approximate mirror symmetry 
about the x axis. Similar "double- well" eigenstates were 
observed for several other oval defects in this sample. All 
of them exhibit the same sequence of states, while the 
number of confined states was varying. 

Localization of the states close to the center of the 
defect, and the absence of mixed states of different par- 
ities indicate a very high potential barrier in the mid- 
dle, and deep wells on both sides, with the correspond- 
ing potential in y direction could be written as V(y) oc 
S(y/b) -fa/ (\y/b\ + 1) with the scaling constants a, b. The 
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FIG. 7: Polariton spectrum of PD3 with Huj c 
detailed description as for Fig.[3] 



1.4353 eV, 



resulting states of the two sides y ^ do not mix signif- 
icantly. The states for y < are shifted by < 1 meV to 
higher energies. In both wells we observe a ground state, 
followed by the first excited state some 4 meV above hav- 
ing one node and the second excited state some 2meV 
further having two nodes, the third excited state some 
2 meV with 3 nodes, and higher states. 

The structural characterization of PD3 is given in 
Fig. [8] The surface profile is oval, with the extension 
along [110] of 3 /im, reduced by a factor of two compared 
to the one of PD1, while the extension along [110] of 6 /am 
is similar to one of PD1. The height of the surface modu- 
lation is about 140 nm, twice the value seen for PD1 and 
PD2, and about twice the nominal height of the GaAs 
Bragg layer. The height increase is consistent with the 
reduced lateral size when accommodating the same vol- 
ume. The defect source is in the 20th DBR period below 
the cavity. 

We now move to the second defect PD4 on MC2, for 
which the hyperspectral transmission images are given 
in Fig.[9j It shows the deepest localized states of all PDs 
studied, with the ground state 40 meV below the contin- 
uum. The states show an approximate mirror symmetry 
along the x and y axis. We can model them with Math- 
ieu functions as shown in Fig.[9| using a strong ellipticity. 
The structural characterization of PD4 is given in Fig.fTo] 
The surface profile is oval as PD3, but with a 20% smaller 
extension and a three time smaller height. The defect 
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FIG. 8: Structural characterization of PD3, detailed descrip- 
tion as for as Fig.[4] 



source is in second period of the DBR below the cavity, 
and has a height of about 100 nm. Being so close to the 
cavity, the additional GaAs has a strong influence on the 
polariton states, and the crater in the middle gives rise to 
a confinement potential with a barrier between the cen- 
ter and the peripheral area, as evidenced by the spatial 
distribution of the confined wavefunctions. 



C. Potential Reconstruction 

The observed localized polariton states can be related 
to an effective confinement potential V m (r) for the in- 
plane polariton motion. We can estimate ^(r) using the 
spectrally integrated density of states D m (Y) created by 
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FIG. 9: Polariton spectrum of PD4 with hw c — 1.4343 eV, detailed description as for Fig. [3] 



Vm(r) below the continuum edge as introduced in Ref.HO 
We use 



D m (r) = £ |* n (r) 



(3) 



n=l 



where the bound state probability densities |\I/ n (r)| 2 are 
taken as the normalized measured intensities 



l*»(r)| S 



J I(r,uj n )dr 2 ' 



(4) 



This expression assumes that the emitted field is pro- 
portional to the polariton wavefunction, which is valid 
for a cavity lifetime which is constant for the in-plane 
wavevector components of the bound states. This is a 
good approximation for small in-plane wavevectors, less 
than 10% of the light wavevector in the cavity of about 
26/fim. Some of the strongly localized states in our study 
with small feature sizes are likely to deviate from this ap- 
proximation. D m (r) is given by the integral of the free 
density of states from zero kinetic energy at the potential 
floor to the continuum when neglecting the spatial varia- 
tion of the confinement potential, i.e. in the limit of small 
level splitting compared to the confinement potential. In 
two dimensions the density of states is constant and given 
by D 2 r> = m/(2nh 2 ), and we find V m (r) = —D(y)/D 2 v- 
We use the effective mass of the polaritons from the mea- 



sured dispersion m = 2 • 10 



where m e is the free 



electron mass. The resulting confinement potentials for 
the investigated PDs are shown in Fig.[TT] The symme- 
try of the potentials reflect the symmetry of the localized 
states. 



D. Surface Reconstruction 

The series of SEM cross-section images S{x 1 y ni z) 
taken at various y n (see Figs. [4j [6j [8| [l0| , provide volume 
information of the defects. To reconstruct the shape of 
the defects in three dimensions, we determine the posi- 
tion of the interfaces between the GaAs and AlAs layers 
in the SEM image. We use PD3 (see Figj8| as an ex- 
ample here. The SEM images show a signal 5, which is 
proportional to the detected secondary electron current, 
differing by about 10% between AlAs and GaAs surfaces. 
The RMS noise in S was about 5% of the GaAs signal. 
SEM images were taken with a nominal magnification 
between 65600 and 79500. We calibrated the vertical (y) 
axis to match the nominal Bragg period, yielding pixel 
sizes between 10 nm and 12 nm with 2% error. The noise 
of S is limiting the precision with which the layer inter- 
face positions can be determined. To enable a reliable 
fit of the interface positions, we have averaged the data 
over 5 pixels (60 nm) along x, orthogonal to the growth 
direction z. The resulting S(x,y n ,z) was fitted with a 
model function S m (z) of the epitaxial structure along z. 
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FIG. 10: Structural characterization of PD 4, detailed descrip- 
tion as for as Fig. [4] 
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FIG. 11: Potentials V m (r) of PDs calculated using Eq.Q. 
The color scale is given, covering to -24meV for PD1 and 
PD2 and to -46 meV for PD3 and PD4. 
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FIG. 12: Example of a fit (blue line) S m (z) to the SEM profile 
S(z) (black line) of PD3, position (x,y n ) = (1.5 /xm, 2.5 fim) 
in Fig.|8] A linear offset has been subtracted for better visi- 
bility. 



The model assumes a Gaussian resolution of the imaging 
with a variance of r/y/2, and a constant thickness d of 
the A c /4 AlAs layers, not affected by the defect, which is 
motivated by the small surface diffusion length of AlAs 
compared to GaAs. The sequence of L AlAs layers in 
GaAs is then described by 



S m (z) = aQ J ra 1 z J ra 2 z + Ax 

L 



(5) 



i=i 



erf 



z - Zi 



erf 



z — z\ — d 



The polynomial coefficients ao,i,2 describe the back- 
ground, A is half the signal difference between AlAs and 
GaAs, and z\ are the positions of the lower interfaces 
of the AlAs layers. The layer index I is the AlAs layer 
number in growth direction. The topmost 3-4 layers were 
excluded from the fitted region as the background varied 



strongly due to the change in secondary electron collec- 
tion efficiency (see e.g. Fig. [4)3). The resolution param- 
eter r was 40 nm corresponding to a FWHM of 67 nm. 
The fitted layer positions z\ show a noise of a few nm. 
An example of such a fit is given in Fig. 12 Using the 
zi(x,y n ) for the different cross-sections n, we can recon- 
struct height maps of the AlAs layers within the struc- 
ture across x and y. A linear slope and an offset along 
x were subtracted from each zi(x,y n ) of an individual 
cross-section n to reproduce the nominal position out- 
side the defect. 

The height maps of PD3 reconstructed from 15 cross- 
sections (see lines in Fig.[8j not all shown) are displayed in 
Fig. 13 The evolution of the surface modulation can be 



followed. The first layer above the defect source (1 = 7) 
shows the center depression of the GaAs, similar to what 
observed in liquid droplet epitaxy 2 -^. Two maxima of the 
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FIG. 13: Height maps of AlAs layers in PD3. The sequential 
numbers of the layers n are given. Linear grey scale from 
-20 nm (black) and +140 nm (white) relative to the nominal 
layer position. The layer 7 is the first above the GaAs layer 
containing the droplet, layer 27 is the first layer above the 2 A 
GaAs cavity layer, and layer 46 is the last fitted layer. 



thickness are observed along the preferential surface dif- 
fusion direction [110]. With increasing layer number, first 
the depression disappears (/ = 9), followed by a general 
extension and flattening of the structure. By integrat- 
ing height profiles of defect for different cuts we have 
determined the volume of the additional GaAs material 
as constant within the error. From this volume, we can 
deduce the radius of the deposited Ga droplet which is 
86 ± 6 nm. 



V. DISCUSSION 

For all of the 15 PDs investigated in this work, of which 
four have been shown as examples, we find a similar ori- 
gin - a local increase in a GaAs layer thickness with a 
depression in the center. The additional GaAs volume 
can only be created by a local deposition of Ga, as the 
growth is limited by the group III element, while the 
group V element As is provided in a much larger amount 
given by the V/III flux ratio of about 50, and desorbs if 
not bound to the surface with a Ga atom to form GaAs. 
The only available source for this excess Ga deposition 
are Ga droplets from the Ga source. 

The shapes of the polariton potentials created by the 
PD are a consequence of the Ga droplet size, and its de- 
position position relative to the cavity layer. In order 
to simulate the 0-dimensional polariton states quantita- 
tively, a full three-dimensional simulation of the mode 
structure in the cavity would be needed, which is be- 



yond the scope of the present work. For a qualitative 
argument, one can use a first-order perturbation picture. 
The polariton intensity is decaying exponentially into the 
Bragg mirror with a decay length of about 400 nm. For 
PD2, the GaAs layer thickening is 50 nm, but it is sepa- 
rated by 23 DBR periods, about 3 /im, or 8 decay lengths 
from the cavity, where the polariton intensity has de- 
cayed to 0.02%. This results in a small influence to the 
polaritons and a small localization energy of the ground 
state of 5 meV. For this defect we observe a large central 
crater and DBR layer discontinuities, as can be seen on 
the central cut in Fig.[4j This could give rise to a re- 
pulsive central part of the PD2 potential as observed in 
Fig.[II] 

In PD1, the Ga droplet had a similar size as in PD2 
(GaAs thickening 90 nm), but it hit the surface only three 
DBR periods above the cavity layer. Even though the 
induced structural perturbation propagates away from 
the cavity layer, is has a much larger influence, with the 
third excited state at -9meV and an estimated ground 
state confinement energy of 20meV. The significantly 
smaller lateral extension of the defects in MC2 leads to 
a larger height of the perturbation, which in turn results 
in stronger confinement of polariton states with shapes 
as seen for PD3 and PD4. 

The evolution of the defect structure during growth 
can be pictured as follows. After a Ga droplet was de- 
posited on the surface, Ga diffuses over the surface from 
the droplet to the surrounding areas. Due the large V-III 
flux ratio, there is sufficient surplus of As2 impinging onto 
the surface to convert the diffusing Ga into GaAs, lead- 
ing to an additional GaAs thickness which decays with 
the distance from the deposition spot, according to the 
Ga diffusion length. The depression in the center of the 
resulting profile is due to reduced GaAs growth below the 
Ga droplet, which requires the diffusion of As through the 
Ga droplet to the GaAs surface. Once the droplet has 
been consumed, the subsequent GaAs growth generally 
tends to smooth the surface due to the Ga surface diffu- 
sion and the preferential attachment of Ga at monolayer 
steps, which have a density proportional to the surface 
gradient for gradients superseding the gradient due to 
monolayer islands (for a island size of 20 nm a gradient 
of 1%). Al instead has a much shorter diffusion length, 
and therefore the surface profile is essentially conserved 
during the growth of the AlAs layers. For Ga grown on 
(001) oriented substrate at 590° C at a V/III flux ratio 
of 2 and a growth rate of 0.25 /im/h the diffusion length 
was reported^ to be 1 /im and 0.02 /im for Ga and Al, 
respectively. 

The observed PD anisotropy of 1:2 to 1:3 along the 
[110]: [110] directions for MC2 grown at a temperature of 
590°C reduces to less than 1:1.1 for MCI grown at 715°C. 
This finding can be explained by temperature depen- 
dent diffusion lengths D for these two crystallographic 
directions. In Ref. 27 £>[no] = 4D[ 110 ] was found result- 
ing in diffusion lengths of L[no] = 2L[ 110 ] for a V/III 
flux ratio of 1.5 and growth temperatures in the range 
of 600° C, in agreement with the aspect ratio of the de- 
fects found in MC2. The reduction of the anisotropy for 
the higher growth temperature of MCI indicates an ac- 
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tivated diffusion with different activation energies in the 
two directions. At higher temperatures, the thermal en- 
ergy supersedes the activation energies and a kinetically 
limited isotropic diffusion is recovered. The presence of 
different activation energies for diffusion in the two crys- 
tallographic directions is plausible as during growth the 
GaAs surface shows a reconstruction^ giving rise to a 
channel-like structure along [110]. 



VI. CONCLUSIONS 

We have shown that oval or round defects in MBE 
grown GaAs microcavities create zero-dimensional po- 
lariton states of narrow linewidths. We have revealed 
their three-dimensional structure and their formation 
mechanism, an impinging Ga droplet during growth. 
While we have deduced effective confinement potentials 
for the defects, a quantitative modeling of the polariton 
spectra from the three-dimensional structural informa- 
tion obtained by the FIB/SEM data is presently missing. 
In the context of polar itonic devices pi our work indicates 
an approach to manufacture two-dimensional polaritonic 
traps by intentional creation of Ga droplets at a spe- 
cific position during the MBE growth of a microcavity, 
rather than ex-situ etching as described in Ref.[30l One 
could also use Ga droplet epitaxy 31 with a low density 
to create well-defined localized polariton states in micro- 
cavities. The narrow linewidths of the polariton states 
formed in this way are favorable for zero-dimensional po- 
lariton switches.^ 
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VIII. APPENDIX 

A. Quantitative Differential Interference Contrast 
Microscopy 

Differential interference contrast microscopy (DIC), 
also know as Nomarski microscopy, was used in reflection 
in this experiment. A Nomarski prism assembly (DIC 
Slider U-DICT with Polarizer U-ANT) is mounted in a 
Olympus BX-50 upright microscope. The illumination 
from a mercury arc lamp is split by the Normarki prism 
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FIG. 14: DIC contrast Cdic (black line) as function of the 
sample position s along the shear direction, and resulting 
height profile h (red line) calculated using Eq.(|8|. 

into two beams 1,2 shifted by the shear displacement A 
in the object plane, with linear polarizations along and 
orthogonal to A. The reflected beams are recombined 
by the prism, creating a polarization state depending on 
their relative phase cpi — (f2- The transmission through 
the polarizer depends on the polarization state, such that 
the intensity depends on the relative phase in the way 

2/dic = I (1 - cos (0o + <fi- ^2)) (6) 

where the offset phase <p is introduced by an adjustable 
spatial offset of the Normarski prism along the optical 
axis from its nominal position for which the beams are 
not displaced in the directional space (objective back fo- 
cal plane). The shear A is similar to the optical res- 
olution of the microscope objective, which allows to ap- 
proximate the phase difference between the two beams in 
first order as the shear times the phase gradient at the ob- 
served position, (pi — (f2 ~ A-V(/?, such that Eq.Q can be 
written as 2/dic = I (1 — cos (<&> + A • Vy?)). Choosing 
<j) Q = ±7r/2, and developing up to first order in the phase 
difference, we get 2/^ IC = J(1±A-Vy?). Measuring 
Jdic f° r both offset phases, we determine the contrast 

C DIC = ffi c ~^ IC =A-Vy (7) 

^DIC + ^DIC 

We can now integrate the contrast along the shear direc- 
tion to retrieve the phase cp. In reflection, the phase is 
related to the surface height h by tp = Airh/X with the 
wavelength A of the light, such that we arrive at 

We assumed here that sample is not birefringent and 
that the phase shift of the reflected light is given by 
the height of the sample surface only, neglecting inter- 
nal interfaces. The latter is justified as the green light is 
absorbed strongly by the structure. The height h(s) was 
determined using Eq.(|8| with s along the direction of the 
shear A. In the measurements presented in this work we 
used a UplanFL 20x/0.5NA objective, for which the shear 
was determined to be |A| = 0.55 /im using a calibration 
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slide in transmission DIC consisting of a PMMA pattern 
of a 200 nm thickness on a glass cover slip, in which case 
the phase is given by <p = 2irh(n\ — 712) with the refrac- 
tive index difference n\ — rt2 = 0.48 between PMMA and 
air. To compensate for systematic errors, the measured 



Cdic across the center of the defect was corrected by the 
Cdic along a line displaced perpendicular to the shear, 
just outside of the defect. An examples of a measured 
Cdic and the resulting height profile h(s) for defect PD4 



are shown in Fig. 14 
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